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Abstract. 

Wc discuss the magnon spectrum of an antiferromagnet (AF) in a magnetic field h 
consisting of two interpenetrating AF sublattices coupled by the exchange interaction 
at T = 0. One-ion easy-plane anisotropy is also taken into account. We calculate 
using the 1/5' expansion the gap in the spectrum which is a manifestation of the 
order-by-disorder effect in this system and the optical magnon mode splitting. Both 
of these phenomena originate from the inter-sublattice interaction. We calculate also 
the gap value at /i ~ /ic in the leading order of the small parameter {he — h) /he using 
the magnon Bose-Einstein condensation treatment, where he is the saturation field. 
By comparing results obtained within these two approaches we conclude that the 1/S 
expansion gives a qualitatively correct result at h he even at large one-ion anisotropy 
but it overestimates the gap value. The application is discussed of these results to the 
actively studied AF of the considered type NiCl2-4SC(NH2)2 (DTN). 
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Figure 1. (Color online.) An antiferromagnet in magnetic field is presented 
containing two interpenetrating antiferromagnetic sublattices shown in different colors. 

1. Introduction 

Quantum criticality of magnetic systems driven by easily controllable external 
parameters like, for example, magnetic field and pressure, and properties of frustrated 
magnets have attracted much attention in recent years. One of the most remarkable 
properties of frustrated systems with a continuously degenerate ground state is that 
fluctuations in them can make favorable some particular states via the so-called 
'order-by-disorder' mechanism [H [21 El II]- An antiferromagnet (AF) containing two 
interpenetrating AF sublattices shown in figure [T] is an example of a system of this 
kind if atoms from one sublattice are located in a zero molecular field of atoms 
from another sublattice. It is the case, for example, if sublattices are coupled by 
exchange interaction. It was shown first by Shender [1] that quantum fluctuations 
lead to anisotropic corrections to the energy and collinear orientation of AF sublattices 
is selected in this way via the 'order-by-disorder' mechanism. This anisotropy is 
accompanied by a gap in one of the two initially degenerate Goldstone magnon modes. 
This gap was observed experimentally in the garnet Fe2Ca3(Ge04)3 [5]. 

Quantum critical properties of another magnet of this kind, NiCl2-4SC(NH2)2 
known as DTN, have attracted much experimental and theoretical attention recently 
illTllHlEllinillllinilTSlITllISllTe]. DTN is characterized by the /4 space group 
with a body-centered tetragonal lattice that can be represented as two interpenetrating 
tetragonal AF sublattices and whose Hamiltonian has the form 

^ = "H-SLi + 'HsL2 + ^ini, (1) 



i,j i i 

where T-LsLi is the Hamiltonian of the l-th AF sublattice, I = 1,2 and Si^i is the spin 
operator at site i in l-th sublattice. The inter-sublattice interaction 'Hint is considered to 
be small and is ignored in the majority of works devoted to DTN. Parameters leading to 
satisfactory agreement between theoretical calculations using Hamiltonian ([T]), ^ and 
the majority of experimental data are S = 1, D = 8.9 K, = 2.2 K and J^^y = 0.18 K 
^llj . where subscripts x, y, z denote exchange constants along corresponding axes shown 
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in figure [T} Due to the strong one-ion easy-plane anisotropy the ground state at h = 
is disordered and all spins are predominantly in the = state. That is why there 
are two quantum critical points in the magnetic field applied parallel to the hard z-axis. 
The first one which is located at h = he = 2.1 T separates a disordered and a canted 
antiferromagnetic phases. The second quantum critical point is at h = he = 12.6 T, 
where he is the saturation field. It is between the canted AF and a coUinear ferromagnetic 
phases. 

A recent ESR experiment at < h < 0.8he demonstrates the need to take into 
account a small inter-sublattice interaction in DTN [6]. The model Q, ^ with T-Lint = 
has a doubly degenerate gapless spectrum with one acoustic and one optical mode. It 
was found experimentally in [6] that one of the two degenerate Goldstone modes has a 
gap and the optical mode is slightly split. To illustrate the idea that the inter-sublattice 
interaction is needed to describe the full set of experimental data, the authors of [6] 
introduced a Dzyaloshinsky-Moria (DM) interaction between sublattices and described 
the gap qualitatively by a mean-field theory. The authors pointed out that an exchange 
interaction between sublattices should also be taken into account but they ignored it 
because it complicates calculations considerably and does not lead to the gap on the 
mean-field level. Thus, it looks desirable to reconsider the spectrum renormalization 
in the AF canted phase by taking into account the inter-sublattice exchange coupling 
because it also leads to the gap via the 'order-by-disorder' mechanism, as is described 
above, and the exchange coupling is normally much larger than the relativistic DM 
interaction [T7] . 

This reconsideration is one of the present paper's aims. We discuss below spectrum 
renormalization of the model ([T]), (|2]) at T = in the first order of l/S" expansion and 
by using the magnon Bose-Einstein condensation (BEG) theory (at h ^ he) with the 
following inter-sublattice interaction: 

T^int = yi.jSi,lSj,2 (3) 

which is considered to be small, V J + D. Using the 1/S expansion we obtain 
the optical mode splitting and expression for the gap which is a generalization of the 
well-known formula by Shender derived for h = [4J. The optical mode splitting is of 
zeroth order in 1/S whereas the gap is of first order in 1/S, being the result of quantum 
fiuctuations. These results are inapplicable at small h in the case of large D, when the 
ground state does not have the Neel order. The gap value is derived also at h ^ he in 
the leading order of the small parameter [he — h)/he using the magnon BEG treatment. 
We demonstrate by comparing results obtained within these two approaches that the 
1/5 expansion reproduces qualitatively the field dependence of the gap at h he even 
at large one-ion anisotropy D but it overestimates the gap value. Particular estimations 
using the results of the 1/S expansion show that V ~ 0.1 K is needed to describe the 
optical mode splitting and the gap obtained in DTN in recent ESR experiment [6] at 
h ^ O.She. 
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This paper is structured as follows. We analyze in section |2] the Hamiltonian ([T])- 
(|3| using the 1/S expansion. The neighborhood of the point h = he is discussed using 
the magnon BEC theory in section [3j Application to DTN of the results obtained in 
section [2] and section [3] is discussed in section HI Section |5] contains our conclusions. 
Two appendices are added with some details of calculations. 



2. Spin- wave analysis 

We assume in this subsection that the exchange coupling of AF sublattices is small: 
V D + J. In this case spins in AF sublattices order parallel to each other in zero field 
|lj and they cant opposite to the field direction by an angle 6 at finite h smaller than 
its saturation value he as it is shown in figure [TJ It is convenient to introduce a local 
coordinate frame {x', y', z') in each lattice site in which the mean spin value is parallel 
to 2;'-axis. Spins components in the laboratory coordinate system (x, y, z) are expressed 
as follows via those in the local coordinate frames: 

ox ox' 

Sn,i = Sn,i COS 6 + S^^^ exp(ikoR„,i) sin 6, 
Sn,i = - Sn,i exp(ikoR„,i) sin 9 + S^'^ cos 9, 

Sn,2 = Sn,2 COS 9 + cxp [iko(R„,2 " Q)] siu 9, 

Sn,2 = - Si,2 cxp [iko(R„,2 " Q)] slu 9 + S^^ COS 9, 
where Q = (1/2,1/2,1/2) is a vector connecting a spin of one AF sublattice to the 
neighboring one of another AF sublattice, ko = (vr, vr, vr) is AF vector of an AF sublattice 
and we put the distance between neighbor spins in an AF sublattice to be equal to 
unity. Imaginary exponents in ^ describe the Neel ordering in the XY-plane. We use 
the Holstein-Primakoff representation of spins components in the local coordinate frame 
written as 



Sn,i + iSii = V2SaiJl - alaj2S ^ (l - 4 a„/45) 



Sn,i - ^S^n,i = ^ " «««n/2'5 " a„ ^ ^25 (l - 4 a„/45) 

Sn,i = — S + al^an, 

3^2 + ^Si, ^ V2Sbl (1 - blbjAS) , (5) 
<2 - ^ (1 - b%/4S) b^, 

^n,2 = — S + bl^bn, 

where operators a and b are introduced for the first and the second AF sublattices, 
respectively. Substituting Q and ([s]) into the Hamiltonian ([T|-([3]) and producing the 
Fourier transformation we obtain 

6 

n = J2'Hi, (6) 

i=0 
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where 
1 



N' 



no = SU (2Jo + 2D + K))cos'^ 



2^cos^ 



2Jr 



5(2 Jo + 2D + Vq) cos9-h sin 9, 
U2 = 5]]4ak [^kcos^^ + S{Jo + D){1- Scos^e) - SVoCos'^e + hcose 

k 

s 



(7) 
(8) 



+ - 



s 



+ 9 SZ^^k^k + bl^a-k) [Vk(l + cos^^) + llc-kosin^ 



(9) 



+ - 



S 



J2 («-k^k + a^felk) (^k - Vk-ko) sin^e, 



N is the number of spins in an AF sublattice, 

fi'k = S{Jo + Jk), 
D 



D (1-1/23), 
Ok = ^= ^ fln exp(zkR„,i' 



(10) 

(11) 



&k = ^ 5^ &n exp [zk(R„,2 + Q)] , (12) 



Jk = ^ Jr exp(zkr), Vy, = ^ K' exp [ik{v' + Q)] 



(13) 



vectors r connect sites inside one AF sublattice and vectors r' connect a site of one 
AF sublattice with neighboring sites of another AF sublattice. Terms I-L3 and 7^4 are 
also important for our consideration. We present them in [Appendix A due to their 
cumbersomeness. One obtains the following expression for Vk in the particular case of 
exchange couplings Vu and Vd with four upper and four lower (relative to the z axis) 
neighboring spins: 



46**^^ cos — COS — 

2 2 



[V^ + Vd) cos^ + z(K 



Vj) sin 



(14) 



In particular, it is seen from (14) that Vko = 0. 

It should be pointed out that the one-ion anisotropy term in (|2| is a constant when 
= 1/2 so that all contributions to observable quantities from the one-ion term should 
vanish if 5 = 1/2. In order results obtained within each order in 1/5 to be in accordance 
with this requirement one has to regroup the 1/5* series by taking into account some 
higher order terms in each order in 1/5*. Such a regrouping of the 1/S series can be 
done quite easily in our case. Simple analysis shows that along with each term in 7/0,1,2 
proportional to D there is the same term multiplied by —1/25". As a result the quantity 
D given by (11), which is zero if S" = 1/2, appears in ([7)-([9| instead of D. [TH] 

Dependence of the canting angle on the field in the zeroth order oi 1/S can be 
found in two equivalent ways: by minimization of the classical ground state energy T-Lq 
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cos I 



h/hc, 
1, 



or by putting to zero linear terms Tii. As a result we have 

if h < he, 
if h > he, 

he = 2SJo + 2SD + SVo, 

where he is the classical saturation field. 

It is convenient to introduce the following Green's functions: 

Gaik) = - (ttkal), Fa{k) = -(a^fc4), 
Gt{k) = - {hbl), F,{k) = -{bl,bl), 
G^'ik) = - (6,4), F^{k) = -(6^,4), 

where k = (a;,k) and a, is the Fourier transform of ak('r 
Dyson equations for these Green's functions: 

Gaik) = Goik) [I + ^kGaik) + UMk) + S^G^(fc) + n^F^(A;)] 
Fa{k) = Goi-k) [y^Mk) + U_-,Ga{k) + S^F^(A;) + n^G^(A;)] 
G^(fc) = Go{k) [^"^-kGaik) + U'^kFaik) + + UkF^'ik)] , 

F'^ik) = Go{-k) [S Va(A:) + Il'l-kGaik) + S_fcF^(A;) + n_^,G^(A;)] , 

where GQ{k) = l/{uj + i6) is the bare Green's function, S, and 11, 11^ are normal and 
anomalous self-energy parts, respectively, and k = (— ci;,k). One has a set of equations 
similar to (20) for Green's functions Gfy{k), Fi,{k), G^{k) and F^{k) which we do not 
present here. Expressing h via 6 using (15), we obtain from the bilinear part of the 
Hamiltonian ^ for the self-energy parts in the zeroth order in 1/S 

Sofc = 9k cos^O + {SJo + SD) siia.^9, 
Hnt = sin^e - {SJo + SD) sin^O, 



(15) 
(16) 



(17) 
(18) 
(19) 

We have the following set of 



(20) 



l-Ofc 



(21) 



n 



Ok 



[SVi^il + cosH) + SV^^k,sin^ 
{SV^-SV^-^,)sin''e. 



Notice that these bare self-energy parts do not depend on frequency. General solution 
of (20) is very cumbersome. In the meantime, it can be simplified considerably when 
the following equalities satisfy: 







— Sfc 








= Tfk 




sr. 




— ^k 


— ^k 











(22) 



which hold exactly in two important for our consideration cases. First, for bare self- 
energy parts (21). Second, at A; = (0, ko) and k = (0, 0). It is shown below that k = ko 
is the point at which the gap opens due to interaction between AF sublattices. As our 
aim is to calculate this gap in the first order in 1/S, we assume below that (22) hold. 
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One obtains as a result of straightforward solution of (20) using (22) and the smallness 
ofV 

^l + nl) + OiV')]/V{k), (23) 

(24) 



Ga{k) = [{uj + Ek){uj^ 

F^{k) = [-u.iu' - j:l + ul) + 0{V')] /V{k), 



G^{k) 
F^(fc) 



[{u + Sfc)2 + n^] - 2Uk{uj + s,)n, + o{v')j /v{k), (25) 

-nliio' -Tl- III) - 2S,n,S^ + 0{V')\ /V{k), (26) 



where 



V{k) = {u'-T.lf + Yil + \{Yi 



y2 



V\2 
k ) 

in 



V\2 



- IS 



V\2\ 



+ 4n,E,(n;sr + nrE^). 



(27) 



We have from (|21|) and (|27|) for the bare spectrum square 

2 sr^2 



-±,k 



nL + ls 



V \2 
Ok I 



n 



V \2 
Ok I 



± 



4|SokS^k 



nokn^i^p — isj^^n^* 



^Ok-^-^Ok 



yV*Y[V 12 
^Ok-'-'-Okl 



.(28) 

Neglecting interaction between AF sublattices one has in (28) S^,^ = 11^,^ = and we 
lead to the doubly degenerate (due to two equivalent AF sublattices) spectrum of AF 
in magnetic field having the well-known form 



Ck = S^ (Jo + Jk) Uo + ^k COS 2^ + 2Dsin^^ . 



(29) 



The minimum of this spectrum is at k = ko and the spectrum is gapless at that point 
because Jkg = —Jq. The inter-sublattice interaction splits these two degenerate branches 
in the classical spectrum. In particular, it leads to the splitting of two branches at k = 
which can be observed, in particular, in ESR experiment. The value of this splitting 
can be found from (21) and (28) with the result 

2Jocos2^ + 5sin2^ 



6e 



k=0 



e+,0 - e-,0 



(30) 



'Jo^cos2^ + Jo^sin^^ 

However, the inter-sublattice interaction remains the classical spectrum gapless. 
Really, using (27) we obtain for the spectrum square at k = kg 

i^ko±KY (31) 
-Iloko and Sjko = -n^ko so that e±,ko = 0. 



V \2 

ko) 



and it is seen from (21 ) and (31 ) that Soko 
We show now that the spin-wave interaction leads to the gap in one of these branches. 

Diagrams leading to the gap are shown in figure |2j Corrections from them can be 
expressed as follows: 



kg 



k 



^Ok 



nok)'e 



k ' 



(32) 



where is given by (29). It is seen from (32) that the first nonvanishing correction at 
k = ko is of the order of V"^. One concludes from (21), (31) and (32) that there is one 
Goldstone mode and the mode with the gap A for which we have 



= sin^ e{Jo + D)^Y1 l^kr(Sok - Hok) 



(33) 
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Figure 2. Diagrams contributing to the gap in the first order in 1/S. SoUd and 
dashed hnes denote Green's functions Ga{k), Fa{k) and Gb{k), Fb{k), respectively, 
given by (f7l and (f8). Composite hnes containing sohd and dashed parts denote 
G^(fc) and (k) given by (f9). (a) Diagrams leading to the first 1/S corrections to 
the linear term in the Hamiltonian ([s]). They lead to renormalization of the angle 9 that 
in turn leads to renormalization of the bare self-energy parts (21). (b) Corrections to 



self-energy parts from three-magnon terms (A.l). (c) Corrections to self-energy parts 



from four-magnon terms (A. 2 I 



0.2 



0.4 



0.6 





1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 
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— ■ D=J 






D = 2J 






1 


II 1 \ 



> 



0.2 



0.4 0.6 



Figure 3. (Color online.) The field dependence of the gap A{h) induced by the 



exchange interaction V between AF sublattices given by Eq. ( 33 ) 



In particular, we recover from (33) at h = D = the well-known Shender's expression 



. Notice that the gap dependence on magnetic field is not trivial because the summand 



in (33) is a function of h. We plot in figure ^ the value A{h)/V for some particular 



values of D assuming that there is exchange coupling J between only neighboring spins 



in an AF sublattice and Vu = Vd = V <^ J in (lA). 

It should be stressed one more time that results obtained above are applicable for 
D < J so that the system has long-range Neel order a.t h = 0. If D is sufficiently large 
the Neel order is destroyed by quantum fiuctuations at h < he < he and h = he is a 
new quantum critical point. This situation is realized in DTN in which D ^ J. This 
statement can be demonstrated by the following expression for the mean value of the 
ion magnetization found in the first order in l/S* at D ^ J by neglecting V: 



(S') = ^ - ^ E \/ T^\A - ^- (34) 
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For exchange interaction J between only nearest neighbors we have from (34) (5*^') ~ 

5*— 0.16 {D / J){1 — h/hc). It is seen from this expression that he is finite at D > 39 J S^. 

We derive in the next section an expression for the gap at h ^ h^ that is apphcable 
for arbitrary D. It is also shown there that (33) works qualitatively at h ^ he even at 
> J. 



3. Neighborhood of the quantum critical point h = h^ 



At h ^ he our results obtained above using the 1/5 expansion can be supplemented by 
calculations using the magnon BEC theory [121 [20] which allows to obtain observable 
quantities as series in terms of powers of the small parameter {he — h)/he. As it is 
shown in Refs. [121 120], the 1/S expansion works badly ioi S ^ 1 at h ^ he due to 
strong quantum fluctuations which are properly taken into account in the magnon BEC 
treatment. 

Following Refs. [191 EO] we start with the case of h > he when all spins are parallel 
to the z axis. Using Holstein-Primakoff representation one writes [20] 

Sn,i ~ ^2^4(1 - Talan), (35) 

Sn,i = -S + a\an, (36) 

7 =\- Vl - 1/25 (37) 

for the first AF sublattice and similar equalities for the second AF sublattice (with h 
operators instead of a ones). After these transformations Hamiltonian ([T])-([3]) has the 
form ([6]), where "Hi = 0, 

^2 = X] (akflk + ^k^k) (fi-k - /u) + ( aj^fek + fe^k'^-k) 'S'Vk , 



(3J 



a\a\ 



1,2,3,4 



-Ji-2-J'S{Ji + Ji) + D 



1,2,3,4 



+ a\b\b2b^ Vi + {b\a\a2ai + b\b%a^ V*\m) 



1,2,3,4 



and fi = he — h plays the role of the chemical potential [191 [20]. One obtains in this case 
for Green's functions Ga{k) and G^{k) defined in (17) and (19) 

Ga{k) = {uj-^k)/V{k), 



G^{k) 
V{k) 



/V{k), 



I) 



(40) 
(41) 
(42) 



(w - Sfc - I) [u - Sfe 

Bare self-energy parts in these equations have the following form at /i < (i.e., at 
h > he): 

Sfc = gk- /i, (43) 

SV^. (44) 
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a) A(p) = 



b) p = 



c) Y 



d) 







Figure 4. Vertexes are shown schematically which appear in expression (47) for the 



energy. The meaning of lines (solid, dashed and composite solid-dashed) is the same 
as in figure [2] Plane a) represents an equation for A(p) in the ladder approximation. 



Using (42), (43) and (44) we obtain the magnon spectrum containing two branches 

/cq, where k = k — ko. 



e±,k = S-k 

which is stable because oc and Vk oc K^KyKz at k 



(45) 



As soon as /i becomes positive, the spectrum (45) becomes unstable at k = ko and 



AF long range order arises. In the magnon BEC treatment this transition is described 
in terms of "condensation" of magnons at states characterized by the momentum ko, 
i.e., (ctko) s-iid (&ko) become nonzero and one has to make shifts 



ako ^ flko + e'^" ^/Npa, 6ko ^ K + e'*" ^/Np^,, (46) 

where pa^h are "condensates" densities and 0a,6 are phases. They should be chosen so as 
to minimize the energy, which has the form at h ^ 



E 



-Kpa + Pb) + ^(Pa + pI) + PpaPb + 27p„p6 COs[2(0„ - 0^)], 



(47) 



where A = 2A(/co), and 7 are vertexes schematically shown in figure |4| Vertexes (3 
and 7 describe interaction between sublattices so that /3 = 7 = 0atV" = 0. Vertex 
A is of the zeroth order in V . It can be found in the ladder approximation by solving 
the equation schematically shown in figure |4|^a), as it was done for AF with one AF 
sublattice [121 ISO]- In particular, (3 can be approximated by its bare value, [3 = Vq, 
which is found from the second term in (39) by putting ki 2,3,4 = ko. Other corrections 



to (3 can be omitted because they are of higher orders in V. As is seen from (39), the 



bare value of 7 is equal to zero. As it is shown in figure |4|^c), 7 is constructed from 
other vertexes. We demonstrate below that 7 is the vertex that gives rise to the gap in 
the spectrum. It can be shown that 7 is negative and of the order of V'^. Some more 
details on calculation of A and 7 can be found in Appendix B 
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+ 




If 



+ 



+ 



Figure 5. Diagrams for self-energy parts at a; = and k = kg of the lowest order in 
/i//ic- Double straight and double dashed lines denote, respectively, condensed particles 
of a and h types and correspond to the last terms in (46). Corresponding vertexes are 
shown in figure |4] 



Minimizing the energy given by (47) with respect to paf, and (\)af, we obtain 



(48) 
(49) 
(50) 



Pa = Pb = P, 

/i = (A + /3 + 27)p. 

Diagrams for self-energy parts at w = and k = ko of the lowest order in p/hc are 
shown in figure [5] Explicitly they can be written as follows: 



-/i + 2Ap + /3p, 
(3p + 47P, 
Ap + 27P, 
Pp. 



(51) 
(52) 
(53) 
(54) 



Substituting these expressions to (31) and using (50) we obtain that one of the spectrum 
branch is gapless and we find for the gap in another branch in the leading order in V 



A 



BEC 



4/i VIWA = - h) \a^MI>^ 



(55) 



As is pointed out above, results obtained in the magnon BEC treatment are more 
rehable than those of the 1/S expansion at ^ /ic if 5 ~ 1, because fluctuations 
are properly taken into account in the BEC theory. It is interesting to compare results 
obtained within these two approaches at /i ~ he. Expression (33) for the gap reproduces 
qualitatively the field dependence of the gap because it also has a linear dependence on 
the field at h h/. 



Ai/5 ~ {he - h) 




nok)2e.3 



(56) 



h=hc 
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But the coefficients in ([55]) and ([56]) are different. We plot the value (A1/5 — 
Abec) /i^i/s + ^BEc) in figure |6]for three spin values, 5" = 1, 3/2 and 2. It is seen from 
figure [6] that results obtained within these two approaches become closer as S rises as it 
must be. Notice also that the 1/5* expansion overestimates the gap value. In particular, 
at S* = 1 and D = we have A1/5 ^ 2Abec- Particular calculations show also that the 
difference between Ai/s and Abec is larger in quasi-low-dimensional systems. 



4. Discussion 



To the best of our knowledge the only compound described by the Hamiltonian ([T])-(|3]) 
which has been studied experimentally in a magnetic field is DTN. Unfortunately our 
results obtained using the 1/S expansion are not entirely applicable to DTN because the 
one-ion anisotropy is very large in this material. As a result the ground state at small 
field is that with = and the Neel ordering arises in the interval he < h < he, where 
he ~ 2.1 T and he ~ 12.6 T is the saturation field. Based on the above consideration of 
the neighborhood of the point h = he one could expect that the results obtained using 
the 1/5* expansion are valid qualitatively near and below he- Unfortunately the data 
of ESR experiment in DTN [B] are available only for h < 10 T ^ 0.8he which makes 
impossible a quantitative description using the expressions obtained above. Thus, we 
can make only a number of qualitative statements. According to the crystal structure 
analysis [21] one could expect Vu = Vd = V in DTN. Equations (14) and (30) with 
~ 0.1 K give the optical mode splitting of the same order as that observed in ESR 
experiment [H] at /i ~ 10 T. Estimations made using (33) show that V ~ 0.1 is required 
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also to describe the gap at ^ 10 T observed experimentally. 



Expression (55) for the gap at h ^ can be useful in further experimental 



discussion of DTN. Values of A and 7 obtained for DTN as it is described in Appendix 
Bjare the following (general expressions for them are very cumbersome): 

^DTN = 4.6 K, 
iDTN 



5. Conclusion 



-O.UK-\ (57) 



To conclude, we discuss using the 1/5* expansion a Heisenberg antiferromagnet 
containing two interpenetrating AF sublattices (see figure [T]) coupled by a small exchange 
interaction V which Hamiltonian is given by ([T|-(|3]). The classical spectrum contains 
a doubly degenerate Goldstone mode if one neglects the interaction between AF 
sublattices. In zeroth order in 1/S the exchange interaction between sublattices gives 
rise to the optical mode splitting at k = 0, which can be measured in ESR experiment 



and which value is given by (30), whereas the classical spectrum remains gapless at 



k = ko. In accordance with the previous result by Shender [1] we observe that quantum 



fluctuations lead to the gap in one of these two modes whose value is given by (33) in 



the first order in 1/S. Equation (33) is a generalization of the Shender's expression for 
a finite magnetic field. 

We calculate also the gap value ai h ^ in the leading order of small parameter 



{he — h)/hc using the magnon BEG theory (see (55)). By comparing this result 
with that obtained within the 1/S expansion we conclude that the latter approach 
gives qualitatively correct result at h he even at large one-ion anisotropy D but it 
overestimates the gap value (see figure [6]) . 



Particular estimations using (30 ) and (33 ) show that V ~ 0.1 K is needed to describe 



the optical mode splitting and the gap obtained in DTN in recent ESR experiment [6] . 



Expression (55) for the gap at h he with values of A and 7 given by (57) can be useful 



in further experimental discussion of DTN. 
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Appendix A. Expressions for "^3 and Ha 

We have for and 7^4 in (|6| 



V NTis = i — ^ sin 6^ cos 6* y^ialala^ — a\a2ai) i 2Jn — 8J1 + lOD — — + Vn 

4\/2 n1 V ^^^^^ ^ 
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(A.l) 



- iJ ^smOcose^ialblb^i + 4«2&i)(V'i - V^i-ko) + h.c, 

1,2,3 

t t /l-2sin2e^ cos'^e,^ 3.2. 
JStiii = ala2'^3a.4 I Ji-i [Ji + ^4) + sm t 



1,2,3,4 



E l + + + / sin^6' Dsiv?6 
(01020304 + O4O3O2O1) ( ^ — Ji + 

1,2,3,4 



4 4 

+ ^ 016^0364 (cos^6'V4_2 + sin^6'V4_2+ko) 

1,2,3,4 

-\y^ (44«3&4 + a\h\h2hi) [(1 + cos^d)VA + sin^^l^-ko] + h.c. 

1,2,3,4 

~5 (44^3^4 + 0J03O2&1) (VI - V"i_ko) sin^^i + /i.e., 



(A.2) 



1,2,3,4 



where momentum conservation laws X]^=i ~ ^^o and X]^=i ki = are implied in sums 



of (A.l) and (A.2), respectively. 



Appendix B. Vertices at h ^ he 

We derive in this appendix vertexes A{uj = 0, p) and 7 diagrams for which are shown in 
figure |4](a) and (c), respectively. Let us start with A(a; = 0,p). The main contribution 
to it is given by ladder diagrams which sum can be found by solving the integral 
equation shown schematically in figure |4]^a). This equation can be solved easily using 
the substitution 



A{u = 0,p) = A + BJp, 



(B.l) 



where A and B are constants. Substituting (B.l) into the equation we obtain the 
following set of algebraic equations on A and B: 



Al x + 



D + TSJo 



By = l, 



(B.2) 



j-A (yil - TS) -TS) + Bil- yil - TS)) = - Q + 



where J-' is given by ([37j) and 

M q Tr^ _j_ 
Jk 



X 



k 



(B.3) 



N ^ 



1 



S Jq -\- S t/k s 



JqX. 



The general solution of the system (B.2) is rather cumbersome. But it can be 



brought into the following compact form in the particular case of 5* = 1: 

;i + Jox) 



A = Jo- 



Jox{3Jo + AD)-2D-Jo 



(B.4) 
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B 



2D + Jqx 



Jnfl 



+ D 



Jox(3 Jo + AD) -2D- Jo 
As a result we have at 5* = 1 from ( |B.l ) and (B.4) 

A = 2A(w = 0, ko) = 2{A - BJq) = Jq- 



2 Jo + 3L> 



'Joa:(3Jo + 4D)-2D- Jo' ^^'^^ 

To calculate 7 we have to take into account diagrams shown in figure |4]^c) which 
involve the vertex A{u = 0, p) obtained above and the vertex denoted by the solid 
circle. Diagrams for the latter vertex are presented in figure |4]^d). The result of their 
summation can be expressed as a multiplication of the bare value of that vertex given 
by the last term in ([39]) by the constant 



1 - 



N ^ 



(B.6) 



As a result we obtain 



7 



N ^ 



32^3 



(^^c^k + 4^A(k))^ 



(B.7) 



If there are different exchange constants J^^ Jy and J^ along different directions 

Jy 7^ Jz) one has to trial the solution of the equation for 



(as in DTN, where Jx ■ 
A(w = 0, p) in the form 

A(u; = 0, p) = A + 5, J,p^ + ByJyp^ + BJ,p^. (B.8) 

The resultant expressions are very cumbersome in this case and we do not present them 
here. They lead to values for A and 7 given by (57) with DTN parameters 5 = 1, 
D = 8.9 K, J, = 2.2 K and Jx,y = 0.18 K 
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